Bistability of Anderson localized states in nonlinear random media 
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We study wave transmission through one-dimensional random nonlinear structures and predict 
a novel effect resulting from an interplay of nonlinearity and disorder. We reveal that, while weak 
nonlinearity does not change the typical exponentially small transmission in the regime of the 
Anderson localization, it affects dramatically the disorder-induced localized states excited inside 
the medium leading to bistable and nonreciprocal resonant transmission. Our numerical modelling 
shows an excellent agreement with theoretical predictions based on the concept of a high-Q resonator 
associated with each localized state. This offers a new way for all-optical light control employing 
statistically-homogeneous random media without regular cavities. 
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The localization of waves in disordered media, also 
known as Anderson localization, is a universal phe- 
nomenon predicted and observed in a variety of classical 
and quantum wave systems Recent renewed interest 
to this phenomenon is driven by a series of experimen- 
tal demonstrations in optics 0] and Bose-Einstein con- 
densates 0]. One of the important issues risen in these 
studies is that the disordered systems can be inherently 
nonlinear, so that an intriguing interplay of nonlinearity 
and disorder could be studied experimentally. 

Nonlinear interaction between the propagating waves 
and disorder can significantly change the interference 
effects, thus fundamentally affecting localization 0, Hj]. 
However, most of the studies of the localization in ran- 
dom nonlinear media deal with the ensemble-averaged 
characteristics of the field, such as the mean field and 
intensity, correlation functions, etc. These quantities de- 
scribe the averaged, typical behavior of the field, but 
they do not contain information about individual lo- 
calized modes (resonances), which exist in the localized 
regime in each realization of the random sample 
These modes are randomly located in both real space and 
frequency domain and are associated with the exponen- 
tial concentration of energy and resonant tunnelling. In 
contrast to regular resonant cavities, the Anderson modes 
occur in a statistically-homogeneous media because of the 
interference of the multiply scattered random fields. Al- 
though the disorder-induced resonances in linear random 
samples have been the subject of studies for decades, the 
resonant properties of nonlinear disordered media have 
not been explored so far. 

In this Letter we study the effect of nonlinearity on 
the Anderson localized states in a one-dimensional ran- 
dom medium. As a result of interplay of nonlinearity and 
disorder, bistability and nonreciprocity appear upon res- 
onant wave tunnelling and excitation of disorder-induced 
localized modes in a manner similar to that for regular 



cavity modes. At the same time, weak nonlinearity has 
practically no effect on the averaged localization back- 
ground. 

First, we consider a stationary problem of the trans- 
mission of a monochromatic wave through a one- 
dimensional random medium with Kerr nonlinearity. 
The problem can be reduced to the equation 



(1) 



where ip is wave field, x is coordinate, k is wave number 
in the vacuum, n — n(x) is the refractive index of the 
medium, and — \ 1S tlr e Kerr coefficient. 

In the linear regime, xl"0| 2 = 0, the multiple scat- 
tering of the wave on the random inhomogeneity n 2 (x) 
brings about Anderson localization. The main signa- 
ture of the localization is an exponential decay of the 
wave intensity, / = \ip\ 2 , deep into the sample and, thus, 
an exponentially small transmission 0, Hiij : ~ 
/in exp(— 2L/l\ oc ) <C 1. Here L is the length of the sam- 
ple and l\ oc is the localization length which is the only 
spatial scale of Anderson localization. Along with the 
typical wave transmission, there is an anomalous, reso- 
nant transmission, which accompanies excitation of the 
Anderson localized states inside the sample and occurs at 
random resonant wave numbers k = fc rcs o [l|-|9j]- In this 
case, the intensity distribution in the sample is charac- 
terized by an exponentially localized high-intensity peak 
inside the sample, / pea k 3> lim and a transmittance much 
higher than the typical one: if^ ^ ^out ( see Fig- !)■ 

Excitation of each localized mode inside the random 
sample can be associated with an effective resonator cav- 
ity located in the area of field localization and character- 
ized by high quality factor Q> 1 [ll| . According to this 
model, the transmittance spectrum T(k) in the vicinity 
of a resonant wavelength, \k — fc res o| ^ fcroso is given by 
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FIG. 1: (Color online.) Excitation of localized modes in the 
linear regime: (a) normalized intensity of the field, I/I m (log 
scale) , vs. the wavelength and position in the sample and (b) 
the transmission spectrum. 



the Lorentzian dependence 0, EH : 



T(k) 



(2) 



+ [2Q (fc/fcrcsO - 1)] 

where T Tes = T(fc rcs o)- 

Let us account now for the effect of weak nonlinearity: 
\x l P 2 \ ^ 1- The nonlinearity becomes noticeable, first of 
all, at the resonant intensity peaks / pca k inside the sam- 
ple. It is physically clear that the Kerr term in Eq. (Q} 
changes the effective refractive index of the medium lead- 
ing to the intensity- dependent shift of the resonant wave 
number: fc rcs0 -> fc r cs(/ P cak)- Since the values of / poak 
and I ut are unambiguously connected, the resonant wave 
number is a function of the output intensity, which in the 
case of weak nonlinearity takes the form: 



^-res(^out) — ^rcsO ~T" 



dk r , 



dl n 



Io 



Substitution Eq. ^ into Eq. © yields 



T(Mo»t) = ^ = - 

An 1 



T 



(3) 



(4) 



[Axiout + sy 

Here we introduced two dimensionless parameters A and 
S, which characterize, respectively, the strength of the 
nonlinear feedback and the detuning from the resonant 
wave number: 

k 



A _2Q d\nk T , 
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dl n 
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FIG. 2: (Color online.) Nonlinear deformations of the trans- 
mission spectra of the resonances 1 and 2 from Fig. I at 
different intensities of the incident wave: numerical simula- 
tions of Eq. (1) (curves) and theoretical Eq. (4) (symbols). 
Light-grey stripes indicate three-valued regions for the high- 
intensity curves, where only two of them (corresponding to 
the lower and upper branches) are stable. 



Equation Q establishes relation between the input 
and output wave intensities, which is cubic with respect 
to / ou t • It has a universal form typical for nonlinear res- 
onators possessing optical bistability [l3|. From Eq. §5§ 
it follows that in the region of parameters: 



A5<0, S 2 >3, \ X \Iin> 



1 



3%/3 \A\T„ 



(6) 



the dependence /out (/in) is of the S-type and the sta- 
tionary transmission spectrum T(k) is a three-valued 
function. In most cases, one of the solutions is unsta- 
ble, whereas the other two form a hysteresis loop in the 
-fout (/in) dependence (see Figs. [2] and [3]). 

It is important to emphasize two features of the equa- 
tions ([I]) and ([5]) , describing the nonlinear resonant trans- 
mission through a localized state. First, they have been 
derived without any approximations apart from the natu- 
ral smallness of the nonlinearity and Lorentzian shape of 
the spectral line. Second, although the resonant trans- 
mission, the effect of nonlinearity, and bistability owe 
their origin to the excitation of the Anderson localized 
mode inside the sample, equations ^ and (0 contain 
only quantities which can be found via outside measure- 
ments. Indeed, T rGS , k lcs0 , and Q are determined from 
the transmission spectrum in the linear regime, Eq. 
while the derivative d In k rcs /dI out | 7 _ can be retrieved 
from the shift of the spectral line when the intensity is 
changed. This enables one to obtain the whole depen- 
dence /out (/in, k) for any given resonance performing ex- 
ternal measurements of T(k) at only two different inten- 
sities of the incident wave. 

To verify theoretical predictions, we numerically model 
the transmission of light incident from x — through a 
random sample consisting of N = 19 alternating layers 
with dielectric constants n\ — 1 and n\ — 10 and random 
widths uniformly distributed in the range (0.12, 1.08)/im, 
Fig. 1. This corresponds to ^ oc = 3.53/xm. Figure [5] 
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shows nonlinear deformations of the resonant transmis- 
sion spectra T(k) for different values of xTn > 0, which 
exhibit transitions to bistability. The analytical depen- 
dence T(k) given by Eqs. ^ and (JS|) with the parameters 
T Tes , Q, and A found from the numerical experiments are 
in excellent agreement with the direct numerical solutions 
of Eq. ([1]). In numerical simulations of stationary regime 
we used the standard 4-th order Runge-Kutta method. 
We note, that the incident field amplitude is a single- 
valued function of the transmitted field. Thus, we solve 
second-order ordinary differential equation Eq. ([T]) using 
transmitted field as the initial condition for the equation. 

The dimensionless parameters T res , and Q from Eqs. 
(jH [5]), can also be estimated from a simple resonator 



model of the Anderson localized states 0, S 11 \- 



AT,T 2 

-L res — , . o 

(Ti + T 2 ) 2 



4^rcs0^1oc 



(7) 



where 



Ti — CXp [-22!resAoc] , T 2 ~ exp [—2{L - Xres)/iloc] 

(8) 

are the transmission coefficients of the two barriers that 
form the effective resonator, a; res is the coordinate of the 
center of the area of field localization, l\ oc is the localiza- 
tion length, and L is the length of the sample. 

Introducing a weak Kerr nonlinearity into the res- 
onator model, one can also estimate the nonlinear feed- 
back parameter A: 

A~Q/T 2 tf, (9) 

where n 2 is the mean value of n 2 (x). 

It is important to note that each disorder- induced res- 
onance is associated with its own effective cavity, so that 
the disordered sample can be considered as a chain of 
randomly located coupled resonators [14] . 

Equations (7)-(9) enable one to estimate the values of 
the parameters describing the nonlinear resonant wave 
tunnelling in Eqs. (4) and (5) by knowing only the basic 
parameters of the localization - the localization coordi- 
nate and the localization length. In particular, substitut- 
ing Eqs. (7)-(9) into Eq. (6) and taking into account that 
the most pronounced transmission peaks correspond to 
the localized states with x ~ L/2 and T\ ~ T 2 , we esti- 
mate the incident power needed for bistability of localized 
states: 



1x14 



> 



exp(-2L/li oc ) 

fcresO^loc 



(10) 



For the parameters used in our simulations this gives 
quite reasonable value |x|T n <) 1CT 5 . If we increase 
the length of the sample, the Q-factors of the resonances 
grow, and the incident power needed to observe bistabil- 
ity becomes smaller. 
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FIG. 3: (Color online.) Stationary and FDTD simulations 
showing hysteresis loops in the output vs. input power de- 
pendence for resonances 1,2, and 3 in Fig.fJJ Panel (d) shows 
deformation of the transmitted Gaussian pulse corresponding 
to the hysteresis switching on the resonance 2. 



To demonstrate temporal dynamics upon the bistable 
resonant tunnelling, we implemented an explicit iterative 
nonlinear finite-difference time-domain (FDTD) scheme 
for modelling pulse propagation through the disordered 
nonlinear sample. For precise modelling of the spectra of 
narrow high-Q resonances, we employed fourth-order ac- 
curate algorithm, both in space and in time 15[, as well 
as the Mur boundary conditions to simulate open bound- 
aries and total- field/scattered-field technique for exciting 
the incident wave. Sufficient accuracy was achieved by 
creating a dense spatial mesh of 300 points per wave- 
length (dx = A/300). To assure stability of the method 
in nonlinear regime, we used the time step of dt — dx/3c, 
and each simulation ran for N — 2 * 10 8 time steps. 
To compare the results of the FDTD simulations with 
the steady-state theory, we consider transmission of long 
Gaussian pulses with central frequencies and amplitudes 
satisfying conditions (fTTTjl . see Fig. Efd). With an appro- 
priate choice of the signal frequencies, we observe hystere- 
sis loops in the Tmt(Tn) dependences which are in excel- 
lent agreement with stationary calculations, as shown in 
Figs.^a-c). Transitional oscillations typical for bistable 
nonlinear structure accompany jumps between two sta- 
ble branches [l6j], and strong reshaping of the transmit- 
ted pulse evidences switching between the two regimes 
of transmission, Fig. EJd). Period of the transitional os- 
cillations is defined by mismatch between external wave 
frequency and nonlinear eigenfrequency, whereas the de- 
cay rate is defined by the Q-factor. We note, that dif- 
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FIG. 4: (Color online.) (a) Non-reciprocal transmission 
through the nonlinear disordered structure, showing differ- 
ent output powers for identical waves incident from different 
directions, (b) Corresponding shape of the incident pulse, and 
pulses transmitted in different directions. 



terministic, and, therefore, these modes can be used for 
efficient control of light similar to regular cavity modes. 
Numerical modelling shows an excellent agreement with 
theoretical analysis based on the concept of a high-Q res- 
onator associated with each localized state. Our results 
demonstrate that, unlike infinite systems, the Anderson 
localization in finite samples is not destroyed by weak 
nonlinearity - instead it exhibits new intriguing features 
typical for resonant nonlinear systems. 

The authors acknowledge a support of the Australian 
Research Council (Linkage International and Discovery 
projects), European Commission (Marie Curie fellow- 
ship), and the Science Foundation Ireland. 



ferent choice of the signal frequencies near the resonance 
can lead to various other behaviours of output vs. input 
curves, with transmission either increasing, when non- 
linear resonance frequency shifts towards the signal fre- 
quency, or decreasing in the opposite case. 

In addition to bistability, the resonant wave tunnelling 
through a nonlinear disordered structure is nonrecipro- 
cical. As is known for regular systems, nonsymmetric 
and nonlinear systems may possess nonreciprocal trans- 
mission properties, resembling the operation of a diode. 
An all-optical diode is a device that allows unidirectional 
propagation of a signal at a given wave leng th, which may 



become useful for many applications [lg, |17|. A disor- 
dered structure is naturally asymmetric in the generic 
case, and one may expect a nonreciprocal resonant trans- 
mission in the nonlinear case. To demonstrate this, we 
modelled propagation of an electromagnetic pulse im- 
pinging the same sample from different sides and monitor 
the transmission characteristics. One case of such nonrc- 
ciprocical resonant transmission is shown in Fig. 4. We 
observe considerably different transmission properties in 
opposite directions with the maximal intensity contrast 
between two directions 7.5:1. Moreover, the threshold of 
bistability is also significantly different for two directions: 
there is a range of incident powers, for which the wave in- 
cident from one side of the sample is bistable, while there 
is no signs of bistability for the incidence from the other 
side. Figure SJb) shows the pulse reshaping for incidence 
from opposite sides of the structure. 

To conclude, we have studied the wave transmis- 
sion through a weakly-nonlinear statistically homoge- 
neous one-dimensional random medium and demon- 
strated novel manifestations of the interplay between 
nonlinearity and disorder. We have shown that even 
weak nonlinearity affect dramatically the resonant trans- 
mission associated with the excitation of the Anderson 
localized states leading to bistability and nonreciprocity. 
Despite random character of the appearance of Ander- 
son modes, their behavior and evolution are rather de- 
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